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Recently, the author has constructed a series of four dimensional non-critical string 
CSJ | theories with eight supercharges, dual to theories of light electric and magnetic 
charges, for which exact formulas for the central charge of the space-time supersym- 
metry algebra as a function of the world-sheet couplings were obtained. The basic idea 
was to generalize the old matrix model approach, replacing the simple matrix integrals 
by the four dimensional matrix path integrals of M = 2 supersymmetric Yang-Mills 
theory, and the Kazakov critical points by the Argyres-Douglas critical points. In the 
present paper, we study qualitatively similar toy path integrals corresponding to the 
^ ■ two dimensional J\f = 2 supersymmetric non-linear a model with target space CP N 
and twisted mass terms. This theory has some very strong similarities with M = 2 
super Yang-Mills, including the presence of critical points in the vicinity of which the 
large N expansion is IR divergent. The model being exactly solvable at large N, we 
can study non-BPS observables and give full proofs that double scaling limits exist 
and correspond to universal continuum limits. A complete characterization of the 
double scaled theories is given. We find evidence for dimensional transmutation of 
the string coupling in some non-critical string theories. We also identify en passant 
some non-BPS particles that become massless at the singularities in addition to the 
usual BPS states. 
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1 Introduction 



In two recent papers Jl|, 0, unexpected properties of the large N limit of Af = 2 
super Yang-Mills theory with gauge group SU(iV) have been discovered and exploited. 
The role of instantons at strong coupling, which has always been elusive, has been 
elucidated in this context by computing the large iV expansion of BPS observables. It 
turns out [flj that the large instantons disintegrate into 'fractional instantons' which 
give non-trivial contributions at each order in 1/N. These 'fractional instantons' are 
thus in particular responsible for the presence of open strings in the string theory 
dual |J, in addition to the familiar closed strings contributing at each order in 1/N 2 
||. The fractional instanton series have a finite radius of convergence, and they 
diverge precisely at the singularities on moduli space. This breakdown of the large N 
expansion is interpreted JTJ as coming from infrared divergences due to the presence 
of a critical point. Similar divergences were encountered long ago in the study of the 
large N limit of ordinary zero dimensional matrix integrals [Q near critical points ||. 
In those simple cases, it was shown in that the divergences can be used to define 
universal double scaling limits from which one can extract exact results in continuum 
string theories. The critical string theories were defined in less than two space-time 
dimensions (the c = 1 barrier) because the only tractable cases were zero or one 
dimensional path integrals. This limitation was overcome in ||, where it was argued 
that the divergences found in [Q] can also be used to define double scaling limits, 
yielding exact results in four dimensional non-critical (or five dimensional critical) 
string theories. The string theories so obtained are dual to theories of light electric 
and magnetic charges which do not have any obvious description in terms of a local 
lagrangian quantum field theory. 

The four dimensional results of [I], were derived by studying the large N ex- 
pansion of the Seiberg-Witten period integrals J7L EL These periods give the central 
charge Z of the supersymmetry algebra as a function of the moduli, and thus the 
masses of BPS states 

M BPS = V2\Z\. (1.1) 

This class of observables is ideal to deduce the strong coupling behaviour of instantons 
because the perturbative quantum corrections stop at one loop and all the non-trivial 
contributions can be understood as coming exclusively from instantons. For our 
purposes, however, the consideration of those special protected amplitudes is not 
enough. To give a full proof of the existence of double scaling limits, one must 
study in principle all the observables, including those with a non-trivial perturbative 
expansion, or equivalently the full path integral. Moreover, the heuristic picture for 
the appearance of a continuum string theory in the limit relies on the observation 
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that very large Feynman graphs dominate near the critical points ||, a fact that 
can in principle be checked on generic amplitudes but obviously not on the BPS 
observables for which perturbation theory is trivial. A related point is to understand 
the universality of the double scaling limits. One argument for universality, put 
forward in 0, is that the double scaling limits are always low energy limits of the 
original field theory. It was observed, however, that to a given CFT in the infrared 
can be associated two different double scaled theories (first class or second class 
singularities in the terminology of Jj]). It would thus clearly be desirable to have a 
more precise characterization of the string theories obtained in the scaling limits. 

The purpose of the present paper is to shed some light on all the above issues 
by studying a particular two dimensional model which is a very close relative to 
M = 2 super Yang-Mills in four dimensions. The model has an exactly calculable 
central charge with the same non-renormalization theorems as in four dimensions 
and the same BPS mass formula (|TTT|) . The analysis of [3], U can thus be reproduced, 
with, as we will demonstrate, qualitatively the same results (appearance of fractional 
instantons, breakdown of the large N expansion at critical points, possibility to define 
double scaling limits for which exactly known BPS amplitudes have a finite limit). 
Moreover, and this is our main point, as our two-dimensional model is solvable in the 
large N limit, we can go far beyond the analysis of the BPS observables. We will 
actually be able to give full proofs of the existence of the double scaling limits, and 
we will give a complete characterization of the double scaled theories. 

The two dimensional theory we consider is the supersymmetric non-linear a model 
with target space CP N and M = 2 preserving mass terms for the would be Goldstone 
bosons. The close relationship of this model with super Yang-Mills was emphasized 
in [fj, and the general analogy between mass terms in non-linear a models and Higgs 



vevs in gauge theories was discussed at lenght in |ID| . It turns out that the two 
dimensional supersymmetric CP N model and the four dimensional super Yang-Mills 
theory are both asymptotically free with a dynamically generated mass scale, have 
instantons, share the same types of non-renormalization theorems, and have both 
BPS solitonic states that can become massless at strong coupling singularities. The 
moduli space of the gauge theory is analogous to the space of mass parameters of 



the non-linear a model |10[ . This very strong analogy can even be made quantitative 



if one adds Nf = N matter hypermultiplets in the fundamental to the pure SU(iV) 
theory, and choose the hypers masses rrii to match the Higgs vevs (pi. One can then 
show || that the central charges of the four dimensional Yang-Mills theory and of 
the two dimensional non-linear a model are actually equal as functions of the m^s, 

Z YM {(pi, rrii = (pi) = Z amodcl (mi). (1.2) 
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The simplest way to understand this relation is to look at the respective brane con- 



structions of the models |Tl|, [12] . The central charges Z are determined by the shape 
of Neveu-Schwarz five-branes which are bent by D4 branes ending on them. It turns 
out that the relevant configurations of branes are the same for the two theories, from 



which ( |1.2| ) follows. In [JT3[| it was argued that equation ( |1.2| ) also implies that the 



stable BPS states are the same in two and four dimensions. It is indeed known in the 
simplest case N — 2 [14 that the knowledge of Z goes a long way toward determining 



the BPS spectra. 

We have organized the paper as follows. In Section 2, we give a rather detailed 
presentation of various classic results about our two dimensional model, including the 
derivation of the central charge as a function of the masses mj. Our goal was to make 
the paper as self-contained as possible. Taking for granted the formulas (|2.28| ) and 



( |2.37| ), the reader may wish to proceed directly to Section 3 where the large N limit of 
the BPS mass formula is analysed, and the double scaling limits are defined. Section 
3 does in two dimensions precisely what was done in [l], |2j in four dimensions, and we 
recover the same qualitative physics (enhangon, fractional instantons, IR divergences, 
double scaled amplitudes). We also discuss at an elementary level the universality 
of the double scaled theories. In Section 4, we give a general analysis of the large 
N expansion. We discuss in details the physical significance of the double scaling 
limits, first in a heuristic way by using the dual Feynman graphs representation, 
then rigorously by using the exact solution of our model at large N. The main 
outcome is a full proof of the existence and universality of the double scaling limits 
exhibited in Section 3. We show that the 'string' coupling undergoes dimensional 
transmutation for first class theories. Another interesting finding is that BPS/anti- 
BPS bound states can become massless at singularities, in addition to the standard 
BPS solitons. We then briefly comment, in Section 5, on other models with M = 1 
or M = supersymmetry, and we conclude in Section 6 by giving possible future 
directions of research. 



2 Classic results 

2.1 Lagrangian, symmetries and renormalization 
2.1.1 J\f = 2 superspace and superfields 

The M = 2 superspace in two dimensions is the dimensional reduction of the standard 
M = 1 superspace in four dimensions, with anticommuting coordinates 9± and 6±, 
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super symmetry generators 



Q ± = Tl — ±2^1, Q± = ±^T29 T d ± , (2.1) 

and supercovariant derivatives 

D± = Tqq- ± 2i6^d ± , D± = ±_ =F 2z0 T <9± . (2.2) 

The two dimensional case has some important peculiarities due to the fact that 
Lorentz transformations do not mix the right and left moving components. One 
can define two R charges, the ordinary fermion number U(1)f under which Q_ and 
Q + have a charge +1, and an axial U(1)a under which Q_ and Q + have respectively 
a charge +1 and —1. Similarly, in addition to ordinary chiral superfields $ defined 
by the equations 

£> + $ = £L$ = 0, (2.3) 
one can define twisted chiral superfields £ by the equations 

D + S = L)_S = 0. (2.4) 

Ordinary and twisted chiral superfields are exchanged by mirror symmetry. Gauge 
fields corresponding to gauge symmetries acting on chiral superfields belong to vector 
superfields V — whose field strengths turn out to be twisted chiral superfields E 
defined by 

E = D + D_V. (2.5) 

The relation ( |2.4j ) and gauge invariance are demonstrated by using the two dimen- 
sional formula 

{D + ,D_} = 0. (2.6) 

In components, the various superfields can be decomposed by using the variables 
y ± = x — 2i8 T 9 T , = x^ =F 2i8 T 8 T , and the Wess-Zumino gauge for V, 

$(x,e,9) = (j)(y) + V26 + ilj--V26-tlj + + 26 + 6^F J (2.7) 
V(x, 6, 9) = -6 + 6 + v + (x) - 0_0_iT + 6 + 6_a + O+ct* + 2i6_6 + {6 + \_ - 0_A+) 

+2z'0+0_(0 + A_ - 0_A+) + 20_0+0+0_£> , (2.8) 
£(z,0,0) = a(y) - 2z(0_A+ + 0+A_) + 20 + 0_(Z> - iv) . (2.9) 



We note that the field strength v = Vqx = DqVi — d±vo is an auxiliary field, a result 
consistent with the fact that gauge fields do not propagate in two dimensions. 
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2.1.2 Lagrangian 



The most general manifestly supersymmetric and renormalizable lagrangian can be 
written as a sum of D-, F- and twisted F-terms, 

L = I Jd 4 8K(<$>,&,E,E j ,V)-Re Jd0-d9 + W d {^)-Re j ' d6„dd + W cl {E) , (2.10) 

where the Kahler potential K is an arbitrary real function and W c \ and W c \ are 
holomorphic functions called the superpotential and twisted superpotential respec- 
tively. The measure of integration over the whole of superspace is defined to be 
d 4 9 = d9 + d9-d9-d9 + . The supersymmetric CP N model is defined in terms of N 
chiral superfields Zi locally parametrizing the complex Kahler manifold CP N with 
Kahler potential 

A N 

^=^ in ( i+ 5>K)- (2 - n) 
y 1=1 

When the coupling g is small, the target space manifold is large and vice-versa. To 
describe CP N globally, we actually need N + 1 coordinate patches 1 < i, j < 
N+l,i^j, related to each other by zf> = z\ k) /zf\ A more elegant description 
of <CP N is in terms of N + 1 complex variables constrained by 

N+l . 

£ w 2 = i < 2 - 12 > 

i=i y 

and with the U(l) identification 



e ia <\> % . (2.13) 

The coordinates z\^ = 4>i/4>j can be used as long as <pj ^ 0. By introducing chiral 
superfields = fa + • • •, a Lagrange multiplier vector superfield V and associated 
£ = D + D_V, and a twisted superpotential 



^. = --rS = -^A + 4) S) ( 2 .14) 
2 2 \2ti g 2 1 ' 



the lagrangian can be written as 

N+l 



1 f f 

L = - / d 4 9 ®U 2V $i -Re d9_d9 + W cX {Y,) . 
•* »=i ^ 



(2.15) 



The vector superfield V implement the gauge symmetry (|2.13|) and the twisted su- 
perpotential implement the constraint ( 2.12 ). The 9 angle term corresponds to the 
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total derivative 9v/(2ir) in the lagrangian. Such a term is actually important at weak 
coupling because of the presence of instantons, for which 

fez. (2-16) 

The 9 term plays a role at strong coupling as well, as we will explain below. By 
integrating out V from (|2.15|) , we recover the pure D-term lagrangian with Kahler 
potential fl2.11| ) plus the topological 9 angle term proportional to 9e tJ,u d i djKd^z i d u Zj. 

As discussed in the introduction, we want to introduce Af = 2 preserving mass 
terms for the would-be Goldstone bosons Zj. These masses play the same qualitative 
role as Higgs vevs in gauge theories. There is no suitable manifestly supersymmetric 
mass term, because a superpotential must be holomorphic and there is no non-trivial 
holomorphic function on a compact complex manifold. However, as first discussed in 
|T5| , non-trivial M = 2 preserving mass terms associated with holomorphic isometries 



of the target space manifold can be written down. An important property of such 
terms is that they induce a non-trivial contribution to the central charge of the 
supersymmetry algebra. In our case, there are N + 1 holomorphic Killing vectors 
associated with the N+1 symmetries 

U(l). : 0.^^%^ ( 2.i7) 

./V of which are independent taking into account the identification (|2.13|) . We thus 
get N independent masses for the N fields Z{. The explicit form of these terms can 
be obtained by gauging the U(l)j symmetries, which amounts to replacing $je 2y $j 
in Q2.15| ) by $je 2( ^ y+vr ^$j, and then by freezing Ej = mj. This procedure of gauging 
also explains why the m^s can be interpreted as the position of branes. To write down 
the form of the final lagrangian, which we will use in Section 4, it is convenient to 
introduce 7 matrices and Dirac spinors, 

7 ° = -a 1 , 7 1 = za 2 , 7 3 = 7 V = a 3 , (2.18) 
V>=(£)» A= (t)' CA= (-A + ) =/i ' ^ = ^°, = ^7°, (2-19) 

in terms of which 

N+1 



l = Y1 (-(^ - ™M (^ + ™")<i>i -w+rrn 
i=i 



2 1 ± 12 



+ipi (ij^dft - 7 M w M + Re(a + m;) - Z7 3 Im(a + m;)) ip, 

Q N+1 4 n N+1 N+1 

+ 2^ v + D^W-^+^fiJ^^-^J^^I 1 - ( 2 ' 2 °) 

i=l " i=l i=l 
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We get the supersymmetric partner of the bosonic constraint ( |2.12| ), J^fJi 1 fifyi = 0> 
as well as the classical potential 



N+l 

Ki = ^k + m,| 2 |0 4 | 2 . (2.21) 

i=l 

The potential yields N+l inequivalent vacua |z), 

4-7J- 

(o-)ixi = -mi , (\4>j\ 2 )i, c \ = —5 , (2.22) 
consistently with the Witten index tr(— 1) F = XEuier (CP N ) =N + 1. 

2.1.3 Symmetries and non-renormalization theorem 

The classical pure CP N model has a SU(iV + 1) x U(1)f x U(1)a bosonic global 
symmetry in addition to M = 2 supersymmetry. The U(1)a symmetry is preserved 
by the twisted superpotential (|2.14|) by assigning U(1)a charge 2 to S. However, 



U(1)a acts chirally on the N + l charged fermions, ipi t ± 1— > e^^ipi^, and will thus be 
anomalous. The anomalous transformation determines exactly the gauge theoretic 
perturbative effective twisted superpotential 

W pcTt (E, m = 0) = ^-+1 S In 4 » (2-23) 
47T eA 

where A is the complexified dynamically generated scale of the theory, 

a n + i = | A |tf+i e tf ( 2.24) 

with a convenient normalization. Similarly, the anomaly can be used to deduce the 
gauge theoretic perturbative twisted superpotential for arbitrary twisted masses m«, 
by assigning charge 2 to the masses, 

1 N+1 S + m 

W> pert (£,m) = — ^(S + m 4 )ln— -p- . (2.25) 

i=l 

These formulas could have been deduced equivalently by a direct computation of 
the quantum corrections to (|2.14 ). This shows in particular that the running of the 



coupling g is given by the a model one-loop contribution, 



1 JV + 1 , w , 

ln rn" ( 2 - 26 ) 



g2(ji) 8tt2 |A| 
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Note that the J\f = 2, N c = Nf = N + 1 gauge theory in four dimensions has precisely 
the same running coupling ||. 

The supersymmetry algebra 

{Q±,Q±} = -4P±, {Q+,Q-} = 2V2Z, {q + ,q_} = o, (2.27) 

implies a BPS bound on the one-particle states masses M > \^2Z, and BPS states 
are defined to saturate this bound fll.l|) . The central charge Z is a linear combination 
of the charges Si associated with the U(l)j transformations ( |2.17| ) and the topological 
charges Tj for solitons interpolating between vacua \i) and \j) (for which T k = 8ik—8jk), 

N+l ^ N+l 

i=i V 2 i=1 

Weff is the exact effective twisted superpotential, whose classical and perturbative 
formulas are given respectively by ( |2.14j ) and ( |2.25| ). The (a)jS satisfy the vacuum 
equation 

dW eS 



da 

and are classically given by (|2.22|) . 



(a= {a) h m) = 0, (2.29) 



2.2 The exact non-perturbative superpotential 

The superpotential ( |2.25| ) has been deduced from an anomaly calculation in the gauge 
theory ( j2.20|) . This gauge theory has an infinite gauge coupling E since there is no 
kinetic term for the gauge field S, whereas the anomaly is computed in perturbation 
theory in E (perturbation theory in E is not to be confused with perturbation theory 
in the non-linear a model coupling g). The formula is nevertheless exact, up to an 
important subtlety discussed at the end of this Section. There are many ways to 



prove this result. One can use the brane construction to compute W e g |12 |. One 



can also use an improved Witten index (T^j to show that Mbps, and thus the central 



charge Z and W e g, do not depend on the D-terms and thus do not depend on E. 
The result (|2.25 ), known to be valid when E — > 0, is thus also valid when E — ► oo. 



The fact that the central charge Z does not depend on E can also be understood 
from Gauss's law, which imply that Z can be computed from the behaviour of the 
fields at large distances, whereas E is an irrelevant coupling in two dimensions. Yet 
another derivation, which is both elementary and rigorous, is to note that the exact 
effective action for S can be deduced by integrating out the fields from ( |2.20| ). 



Since the <3?iS appear only quadratically, this can be done exactly. To isolate the 
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twisted superpotential term from the general non-local effective action for £, one 
uses the fact that the four-momentum P M can be written as an anticommutator of 
the supercharges, and thus that any non-local F-term can also be written as a D- 
term. The most general F-term is then given by the local twisted superpotential. To 
compute W e fi, it is thus enough to consider constant fields, and to set the fermions 
and the field strenght v to zero. The general effective action then admits an expansion 
is powers of D, 

S cS = J d 2 x (-2D Re + 0(D 2 )Y (2.30) 

The .D-terms can contribute only at order D 2 or higher. The term linear in D, 
which is given by a simple one-loop calculation in the gauge theory, together with the 



analyticity properties of the twisted superpotential, yield the formula fl2.25|) . 



There is a difficulty with (|2.25 ), which seems to be at the origin of some confusion 



in the literature. The formula is ambiguous because the logarithm is a multivalued 
function. The ambiguity corresponds to adding a term ipH/2, p G Z, to the twisted 
superpotential, or equivalently to shifting the 9 angle by —2np. If we are in the 
vacuum and if the masses are such that \iTij — >> A for i j£ j, then the physics 
is weakly coupled and the structure of the vacuum is semiclassical. In particular, 
this means that the boundary conditions at infinity are such that the quantization 
condition fl2.16Q holds, and thus the ambiguity on W e s is unphysical. This can actually 
be checked explicitly. The physical content of W e g is summarized by the vacuum 
equation 

dW eB 



da 



0. (2.31) 



Using ( |2.25| ), this reduces to 

N+l 



f[ \a + mj\ = \A\ N+1 (2.32) 



3=1 
N+l 



arg^-^ = 2 P 7T, peZ. (2.33) 

3=1 



When A — > 0, ( |2.32| ) implies unambiguously the classical result ( |2.22j ). For 



\ m 3 



ml 



» A, 



|) and ( j2.33| ) imply a unique convergent instanton series expansion for 



-m, 



{i)f 



)A k(N 



+1) 



(2.34) 



fe=i 
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where the c£'s are calculable functions of the rrijS, for example = 1/ Ylj^i( m j~ m%). 
The important point is that the series ( |2.34j ) does not depend on the ambiguity in 



(pi). 



When the vacuum is no longer weakly coupled, instanton calculus is plagued by 
infrared divergences and the semiclassical approximation is no longer valid. In par- 
ticular, the series ( |2.34| ) has a finite radius of convergence. In the strongly coupled 
regime, the original field variables strongly fluctuates, the classical geometric picture 
of a CP N target space is lost, and the arguments leading to the quantization condi- 
tion ( |2.16| ) do not apply. The (a)iS can nevertheless be calculated, because M = 2 
supersymmetry implies that they are the analytic continuations of the series (21 



The analytic continuations are easily found by noting that those series are the N + 1 
solutions of 



JV+l 



Y[(v + m,)=A N+ \ (2.35) 



and that analyticity implies that fl2.35p is always valid. The unambiguous vacuum 



equation ( |2.35| ) is obtained by integrating 

^ = ^ ln ll^v-' (2 ' 36) 

i=l 

which yields 

i N+1 v , i N+1 v , 

m) = - E In J] — ^ + 7Z E ln ^ ' ( 2 ' 37 ) 



It is interesting to note that this specific formula for W e g corresponds to the lowest 
energy density at fixed a. The qualitative difference between the weakly coupled 
regime where the series ( [2.34 ) converge and the strongly coupled regime where one 



must use the analytic continuations is that in the first instance (a)i(6) = (cr)j(# + 27r 
while in the second instance the different vacua are mixed up when 9 — ► 9 + 2n [[17] 



At strong coupling, we see that the apparent ambiguity in W e g simply corresponds 
to a choice of a particular vacuum. The resolution of the difficulty associated with 
the branch cut in the formulas ( |2.25| ) or ( |2.37| ) is thus qualitatively different at weak 
coupling and at strong coupling (instantons or choice of vacuum), but the physics 
described by ( p. 37] ) is always consistent and unambiguous. 
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3 The BPS mass formula at large N 



In this Section, we start the study of the large N limit of our model, restricting our 
attention to the exactly known central charge ( |2.28| ), in strict parallel to what was 
done in four dimensions in ]T|, [|. In Section 4, we will show that the results obtained 
by studying this restricted class of observables do generalize to the full theory. 



3.1 The enhangon and critical points 

To study the large N limit of the central charge, one must first study the large N 
limit of the roots of the equation (|2.35|) . This problem was solved in Section 3.2 of 



JH, and we briefly recall the results below. It turns out that a consistent large N 
limit is approached when the mass density 



- N+l 

= j^Ti £ 5(2)(m " mi) (3,i: 



»=i 

goes to a well defined distribution p^ when N — > oo. This distribution can be the 
sum of a smooth function and of 5 function terms. Studying the full N+l dimensional 
space of mass parameters, or, at large N, the full space of distributions p^, is not 
very convenient. It is more instructive to consider one-dimensional sections of the 
full space, parametrized by a global complex mass scale v. Given a fixed distribution 
Pn{v) of dimensionless numbers Vi, the masses are defined to be m[ = wi. The 
associated density is then 

(v) f \ 1 ( m \ 
By introducing the dimensionless ratio 



P ( N ) (m) = ^ PN (-). (3.2) 



v 

r = W\ 

and the polynomials 

N+l N+l 



(3.3) 



p(x) = Y[(x + Ui), q(x) = p(x) - e l9 /r N+1 = J[ (x - Xi ) , (3.4) 



i=l i=l 

one can write the central charge (|2.28|) as 

N+l , N+l 



1 

i=l * i=l 



(3.5) 
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where the function 

N+l 



(N + l)x + ^z/ i hi(x + z/ i ) (3.6) 



W{X) = 

1=1 



is the field-dependent part of the twisted superpotential W e s/v ( p.37|) evaluated at 
one of the vacua. 

If Poo is a smooth function, the physics is weakly coupled in all the vacua when 
|r| — ► oo. The "quantum" roots Xi of q and the "classical" roots —Vi of p then 
coincide at large N, up to exponentially suppressed instanton terms. This picture is 
valid as long as \r\ is greater than some critical value r c . When \r\ < r c , the instanton 
series can diverge, and the roots Xi gradually arrange themselves along an inflating 
curve in the x-plane. This curve is a generalization of the enhangon discussed in |fT% 



When r — > 0, the enhangon eventually eats up all the roots, and approaches a circle 
of radius l/|r|. If Poo has 5 function terms, some of the vacua are always strongly 
coupled, whatever large |r| is. The roots corresponding to such vacua are arranged 
on an enhangon for all r. Finally, let us note that the enhangon can have several 
connected components, associated with the connected components of the support of 

Poo- 

Of crucial importance to us are the critical points that are obtained for some 
special values of the mass parameters. These critical points, or singularities, are 



physically similar to the Argyres-Douglas points fL9 |, which were argued in || to 



be the four-dimensional generalizations of the Kazakov critical points found in zero 



dimensional matrix models p0| . Mathematically, both the two dimensional and four 
dimensional critical points are obtained when the discriminant of the polynomial q 
( |3.4| ) vanishes. At large N, it was explained in [||] that this can happen either when a 
classical root is eaten up by the inflating enhangon (first class critical point) or when 
several connected components of the enhangon collide with each other (second class 
critical point). Let us emphasize that the distinction between first class and second 
class does not arise because the low energy physics is different in the two cases — 
the corresponding CFTs are actually the same — but because the large N expansion 
behaves differently near a first class or a second class singularity. In particular, for 
a given CFT, the first class and second class double scaled theories defined in 
are different. In Section 4, we will completely characterize those theories in the two 
dimensional setting of the present paper. 
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3.2 Fractional instantons and IR divergences 

We now give two concrete examples of a first class and a second class singularity in 
our model. The mass densities are chosen to be the same as the Higgs vevs densities 
of the examples studied in [J] . For that reason, some of the formulas of the Section 
4 of [|TJ can be used here. One minor difference is that N in must be replaced by 
N + 1. We have chosen this convention because, in perturbation theory, the iV of 
SU(iV) and of CP N do play the same role, distinguishing between the topology of the 
dual Feynman graphs (see Section 4). 



3.2.1 An example with a first class singularity 

We choose the distribution 

p N (u) = j^j S^(y - 1/(N + 1)) + S&>(v - 1/(N + 1) + 1) (3.7) 

and the 9 angle to be 9 = 7r.[] The first class singularity occurs at the critical 
parameter r c = (iV+ l)/^ 1-1 /^" 1 " 1 ) when the two positive real roots X\ and %2 of the 
polynomial q defined in ( [3.4|) coincide. We want to calculate the central charge of the 
BPS state that becomes massless at the singularity. At large N and |r| > 1, we have 
(see equations (45) and (44) of [|TJ with N replaced by N + 1) 

1 r + ln(r-l) 1 /, . xx2 21n(r-l) \ /n . „. . . 
X ^ = r ]v> + 2iVV ( (ln(r " 1)} " r -1 + 2 V + ° {1/N ) m 

X2 = 1 - l/(A + 1) - l/r N+1 + 0(l/r 2{N+1) ) = 1 - 1/N + 1/N 2 + 0(1/N 3 ). (3.9) 
We then immediately get, by using fl3.5|) and (|3.6|) , 

Amz r-1 , (r- l)(ln(r- 1) - 1) (ln(r - l)) 2 , nl . . 

This formula displays all the important qualitative features of the large iV expansion 
of "instanton generated" BPS observables [jl|: the expansion parameter is 1/N; the 
expansion breaks down at the critical point r = 1; each order in 1/N is given by a 
mixing between In r terms coming from perturbation theory and series in 1/r obtained 
by writing ln(r — 1) = lnr + ln(l — 1/r) = lnr — JZfcLi l/{kr k ). These series in 1/r 
are naturally interpreted as coming from fractional instantons of topological charge 
1/(N + 1). These fractional instantons would be the remnant of the disintegration of 

1 The 9 dependence could be absorbed in the phase of the parameter r. The choice 9 = tt is 
convenient to compare with the formulas of [0 . 
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large instantons at strong coupling M. Let us emphasize, however, that the fractional 
instanton picture remains elusive, because we have not found the corresponding field 
configurations (that must be singular in the original field variables), and also because 
at large N the topological charge is vanishingly small. 



3.2.2 An example with a second class singularity 

Let us now suppose that N + 1 is a multiple of four, choose 9 = 0, and consider 

PnM = ~ (>V " 1) + <5 (2 V + 1)) • (3.H) 

The second class singularity occurs when r = 1, at the merging of the roots 

x x = y/l - 1/r 2 (3.12) 

and X2 = —X\ of the polynomial q. The formula ( |3.12|) is exact. In particular, the 
large N expansion of the roots near a second class singularity, though non-analytic, 
is not blowing up. Unlike the four dimensional case, this implies that the large N 
expansion of the central charge itself is not divergent. The exact formula is easily 
derived, 

= v 1 + Jf) ~ ^ " + ln ( X + V 7 ! " lA 2 )) • (3.13) 

We do not find divergences when r — > 1 because the large N expansion of z has 
only two terms. We will see in Section 4 that the 1/N expansion of more general 
observables has an infinite number of terms and does suffer from IR divergences at 
the critical point. The physical origin of these divergences is actually the same for 
first class or second class singularities. 

It is interesting to note that we get fractional instanton series (presently of topo- 
logical charge 2/(N +1)) in the exact formula ( 3.13Q . The reason is that due to the 



special choice for the density ( |3.11 ), the 1/N expansion has only two terms and gives 



the exact answer in this example. Generically, we expect to get fractional instanton 
series only at large N. This simply comes from the fact that, in order to obtain a 
smooth N —* oo limit, one must take the 1/(N + l)th power of (|2.35|) and write 



exp Jd 2 mp N (m) In (a + m) = A . (3-14) 
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3.3 The double scaling limits 
3.3.1 First class 



Let us consider the scaling 



iV^oo, 5 = r — 1 — > , NS — In iV = constant = 1/k + In k . 



(3.15) 



which was used in similar circumstances in four dimensions (equations (52) and (54) 
of 0). From ( |3.10| ), one can easily show that subtle cancellations make the first three 
terms in the perturbative expansion of the amplitude 



A 



4iirNz 

~7T 



(3.16) 



finite in the double scaling limit ( 3.15Q 



^scaled = 7TT - - + 0(k) . 

1k a k 



(3.17) 



Going beyond the perturbative expansion is actually easy. The exact formula for A 
before the scaling is 



A = N 



dx In 



x + 



N + l 



N / 1 
X + 



N +1 



1 r 



AT+l 



(3.18) 



By changing the variable to u = N — (N + l)x, we immediately see that in the limit 
( p,15| ) we have 



scaled 



e(e-v«/ K ) 



du In • 



D u— 1/k 



1/k 



U U KU 

w in — 

Ik e 



UK 

€(e-V-/K) 



(3.19) 



J 1/k 



where the function ^ is defined by 



v = ue 



u = £(v) for u E [0, 1]. 



(3.20) 



This explicitly demonstrate that the suitably rescaled central charge ( 3.16|) has a finite 
limit in the scaling ( p.!5| ). The function C, gives purely non-perturbative contributions 
proportional to e -1 / K , and thus the perturbative expansion is entirely given by the 
first two terms (|3.17| ). 
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One might wonder to what extent the result ( |3.19| ) is universal, and whether 
one can find generalizations. A straightforward argument for universality, which was 
given in 0, is that the scaling ( p,15| ) corresponds to a low energy limit. Indeed, the 
central charge is related to a mass scale by ( |1.1| ), and equation ( |3.16| ) shows that 
the amplitude having a finite limit in the scaling is N times this mass scale. This 
implies that only the light degrees of freedom, that become massless at the singularity, 
survive in the scaling limit. The result (|3.19 ) suggests that the double scaled theory, 



which must describe the interactions between those light degrees of freedom, is a field 
theory with an effective superpotential 



u 2 u KjU 



w es(u) = — win (3.21) 

2 k e 

We will be able to characterize this theory in Section 4, and universality will then be 
obvious. Right now, we can discuss the dependence of the superpotential ( |3.21 



on 



the particular choice (|3.7| ). It is actually not difficult to treat the general case where 
m classical roots, for example x = . . .x = — u m , melt into the enhangon. The 
starting point is the formula 



N+l 

A = N I dx In 



.N+l 

i=l 



n ( x + ( 3 - 22 ) 



where the upper and lower bounds of the integral are two distinct zeros of the inte- 
grand. The density for the roots on the enhangon 



1 N+l 

d »M = N + 1 - m E ^ (3-23) 

i=m+l 

is taken to be arbitrary as long as it goes to a well-defined distribution d{v) of bounded 
support when iV — > oo. We want to see to what extent the double scaled amplitude 
depends on d(u). The critical point occurs when the classical roots z/j, 1 < % < m, 
melt into the enhangon. As we will see, by adjusting the critical value of r, we can 
assume that the critical value of the z/^s for 1 < i < m is an arbitrary number M 
such that | — M + u\ >e>0 for all v in the support of d[y). Let us thus define 
x = —M + v/N and Ui = M + Vi/N for 1 < i < m. We can write 

N+l m N+l N+l 

Y[[r(x + Vi )\ = + Vi) exp [ln(r N+1 N~ m ) + ^ ln(-M + v t ) + ^ ^ _ M + 

i=l i=l i=m+l i=l 

m 

= J^(u + Ui)exp N(A + lnr) - mlniV + lnr + B + Cwj , (3.24) 

i=i 
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where we have neglected terms that will go to zero when N — > oo, and A, B and C 



are some constants depending on the distribution (|3.23| ) but not on N. In particular, 
the critical value of r is 



exp 



The generalized double scaling limit 



- / d 2 vd{v)\n{-M + v) 



(3.25) 



N -> oo , r -> e~ A , ^ -> — M , JV(z/i + M) = est = = m/C , 

iV(A + lnr) -mln(C'iV) + J B-A-z0 = cst = l//« + mlnK (3.26) 



then yields 

Acaled = g \Weft{u2) ~ W e ff(Ul)J (3.27) 

where Ui and w 2 are roots of the equation 



e 



+ «*) = ( 3 - 28 ) 



1=1 

and 



w eff (it) = — + - + V(u + Wi)to — ^- (3.29) 



1=1 

The parameters itj were chosen without loss of generality such that YlT=i M * = ®- We 
see that all the dependence in the general density ( |3.23| ) we started from is in a trivial 
global finite factor C. In particular, the formula ( |3.21| ) for the case m — 1 is recovered 
after changing u in —u. 

3.3.2 Second class 

The general case of an mth order critical point can be described by choosing iV + 1 
to be a multiple of m and the polynomial p of equation (|3.4|) to be 

p( x ) = [x m + u m , k x m ~ k + 1 J . (3.30) 



fc=2 



Defining 



m— 2 



to = m \n(r N+l e- ld ) , = N l - j/m Uj , T(u) = ^ t k u k + w m , (3.31) 



fc=0 
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and taking the N — > oo limit keeping fixed the tjS, we see that the amplitude 

2mmN l ' m z mN^ m 



A 



2t 



l/m 



dx In 



e~ l \ N+l p{x) 



has a finite limit 



A 



scaled 



2f. 



l/m 



T{u) du , 



(3.32) 



(3.33) 



where the upper and lower bounds of the integral are roots of the polynomial T. This 
formula strongly suggests that the double scaled theory is a simple Landau-Ginzburg 
field theory with a superpotential 



m-2 

E 

fc=0 



t k u 



fc+1 „,m+l 



+ 



U 



k + 1 m + 1 



(3.34) 



For example, if = for & > 1 and to — 1/ K we g e t the roots 

Uj = t^ m exp(i7r(l + 2j)/m) 

and the amplitudes 



A 



K 



,.1/m /•« 



T{u) du = e ^+ k + l )l m sin(7r(j - k) /m) X m (/t) 



with 



1m. ( ^ ) 



[m + 1)k 



(3.35) 



(3.36) 



(3.37) 



Equations ( |3.33j ), ( p.36|) and (|3.37| ) are the exact analogues of equations (41), (42) 
and (43) of 0. 



4 The full large N expansion 

The results of the previous Section suggest that, if we rescale the space-time variables 
from x^ 1 to cr^, 

a » = N -i/P x ^ ( 41 ) 

and take the double scaling limit ( |3.26|) (for p = 1) or (|3.31|) (for p = m), then the 
original non-linear a model tends to a well defined "double scaled" theory describing 
the interactions between the light degrees of freedom. A full proof of this statement of 
course requires to study the full path integral, not only the central charge, and that's 
precisely what we intend to do in this Section. However, before entering into the 
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details, it is useful to give a qualitative discussion that applies to the more difficult 
case of gauge theories as well. 

An important point is that, even though the double scaling limits correspond to 
low energy limits, as (O) clearly shows, the limiting procedure does not introduce 



a cut-off. This means that the resulting theories must be defined on all scales, and 
are thus fully consistent relativistic quantum theories, obtained from an asymptoti- 
cally free quantum field theory by taking a consistent limit. This fact is particularly 
startling in four dimensions, where the double scaled theories are relativistic quan- 
tum theories of light electrically and magnetically charged particles, for which only 
effective descriptions were known. 

A very elegant, if only heuristic, way to elucidate the nature of the double scaled 
theories is to introduce a dual representation for the Feynman graph, and realize that 
very large Feynamn graphs dominate near the critical points. This classic analysis 
[ 21], |2~2"i , that we sketch in the next subsection, suggests that the four dimensional 



double scaled theories are string theories while the two dimensional double scaled 
theories are field theories. We then proceed to an explicit proof of this result in two 
dimensions, where the large Feynman graphs of the original non-linear a model can 
be explicitly summed up. 



4.1 Loops of bubbles and the continuum limit 

A generic observable of the two dimensional CP N model can be expanded at large N 
as a power series in 1/N, 

A = N a Nl ~ h A h , (4.2) 

h>0 

where N a is some normalization that insures that A has a finite limit in the double 
scaling. The coefficients Ah can pick contributions both from Feynman diagrams 
and from non-perturbative effects. In the case of SU(iV) gauge theories, Feynman 
diagrams generate a series in 1/N 2 , while non-perturbative effects can contribute at all 
orders in 1/N [[[]]. In Section 3, we have discussed observables for which perturbation 
theory was trivial. However, many other observables are dominated by the Feynman 
graphs contributions. An example that we will discuss explicitly below is the mass of 
non-BPS states. For those, we can write 

A h = Y,A h , k \ 2k , (4.3) 

where A = g 2 N is the renormalized 't Hooft coupling constant. For the double scaling 
limits to yield a finite result, it is necessary that the coefficients Ah(X) diverge near 
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Figure 1: Part of a typical large Feynman diagram in an SU(iV) gauge theory A 
dual representation (gray lines) is obtained by associating a p-gon to each vertex of 
order p. The p-gons generate a discretized Riemann surface. The double scaling 
limits correspond to a continuum limit where the number of p-gons becomes infinite 
and thus the discretized Riemann surfaces become genuine smooth world sheets. The 
power of k 2 oc 1/N 2 counts the genus of these world sheets. 



the critical points A = A c , at least for sufficiently large h. The whole idea of the 
double scaling limits is actually that those divergences are specific enough so that 
they can be compensated for by taking the N — > oo limit together with the A — > A c 
limit. Typically, one has, up to logarithmic terms, 

A \^ { x-l)^ ' (4 ' 4) 

where 7^ is some susceptibility. This shows that near A = A c , the terms with a high 
power of k dominate in (}4.3|) . Those terms are generically associated with very large 
Feynman diagrams, containing a lot of interaction vertices. 

What do those diagrams look like? In the case of gauge theories (Figure 1), the 
answer |§ is that the diagrams contributing to a given order in 1/N 2 can be mapped 
to discretized Riemann surfaces of a given genus. Large diagrams have a very large 
number of polygons, and thus the double scaling limit is a continuum limit for the 
discretized surfaces. We conclude that the resulting theory must be a string theory. 
For non-linear a models the analogous statements are easy to derive. In the case of 
linear a models [21], the typical large N graphs are "bubble" diagrams, the order in 
1/N being related to the number of loops of bubbles. In a dual representation (Figure 
2), we obtain a discretized world line (or "polymer") with a given number of loops. 
The double scaling limit is then a continuum limit for these discretized loop diagrams, 
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Figure 2: Part of a typical large Feynman diagram in the CP non-linear a model. 
A dual representation (gray lines) is obtained by associating a bound linking p 
"molecules" (small gray disks) to each vertex of order 2p. The bounds generate 
a discretized "polymer" with both forced (for bounds with p > 3) and dynamical 
branching. The double scaling limits correspond to a continuum limit where the 
number of bounds become infinite and thus the discretized polymers become genuine 
smooth world lines. The power of k oc 1/N counts the number of loops (two in our 
example) of these world lines. 

and as a result we should obtain a standard field theory. The case of non-linear a 
models is similar, with the additional subtlety that we have interaction vertices of 
any order p > 2, but we still expect the resulting theory to be a field theory. 

For the purposes of the calculations that follow, it is convenient to go to the 
euclidean for which x 2 = ix° and ^ Eucl = -i^ MW[ . The euclidean lagrangian deduced 
from (P2CD is 

N+l .„ 

L E = -L = J2 ( V ^l + \<r + rrii] 2 \H 2 + & (f ~ o- - m) Vi) + y * v 

i=l 

N+l ^ N+l N+l 

-D (J2 l^l 2 - 4) - E Ml +iV2j2 Mi V , ( 4 -5) 

i=l " i=l i=l 

where we have defined the covariant derivative V a = d a + iv a , the field strenght 
*v = d\V2 — d%Vx, and a = Re a — ia 3 Im a, rftj = Rem; - ia 3 Im m.j. 
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4.2 First class singularities 



We consider the distribution ( |3.7| ) again. It is convenient to make the substitution 
v/(N + 1) and to define 



a 



a 



m 



Ar, 



(4.6) 



in line with the notations of Section 3. The large N limit is studied by integrating 
the superfields $j = (0,,^) from (4J3). This can be done exactly, and yields a non- 
local effective action proportional to N. The 1/N expansion is then a perturbative 
expansion for this non-local effective action. For the particular distribution (|3.7| ), the 
superfield $ n+i plays a special role. One must actually keep this superfield explicitly 
in order to get a well-defined saddle point at large N. We are then left with a path 
integral over the fields = (J)n + i/VN, ip = ^>n+i/VN and E = {a,n,D,v) which 
reads 



V((f>,ip,a,fi,D,v c 



det r (y 



cr 



det r 



-V 2 + \a\ 2 



D -2ji(f -<t)-V) 



N 



cxp 



-N J d 2 x (v a + V Q + (\a + m\ 2 - D) |0| 2 + if, (f - a - rh) ^ 



(4.7) 



The renormalized determinants are studied in Appendix A, where all the formulas 
that we will use in the following can be found. The scale /i appearing in ln(/x/|A|), 
not to be confused with the spinor //, is a renormalization scale appearing in the 
definition of the determinants. 



4.2.1 The BPS/anti-BPS bound state 



Before taking the scaling limit (|3.15 ) on the full path integral fl4.7|) , it is instructive 
to study explicitly an observable that has a non-trivial perturbative expansion, to 
complement the discussion of Section 3. We will consider the mass of a BPS/anti- 
BPS bound state that turns out to become massless at the first class singularity. 

The saddle point equations for (|4.7|) are deduced from the effective potential^ 



V cS = (\a + m\ 2 -D) 



D . n 
— ln-V- + 

2tt A 4tt 



In 



cr 



e/i z 



47T 



D 



2 We have chosen the background electric field to be zero. A possible non-zero electric field is an 
effect of order 1/N, and thus can be neglected in the saddle point equations. 
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The saddle points dV e s = correspond to the possible N = oo limit of the vacuum 
expectation values of the fields, and are also given by the condition V e s = by 
super symmetry. Two cases must be considered. When \r\ = \m/A\ < 1, the only 
solution is 

(0) = (D) = O, |<a)| = |A|. (4.9) 

We thus get an enhangon, as discussed in Section 3.1 or in more details in [I], which 
is simply a circle in the a plane. When \r\ > 1, in addition to the enhangon ( |4.9| ), we 
get another solution, 

! 777 

This solution corresponds to the root x 2 in equation (|3.9| ), while the root x% of equation 
( p.8|) lies on the enhangon (|4.9|) . The critical point occurs when the two roots coincide. 

The vacuum ( [4.1 0|) is weakly coupled when \m\ ^> |A|. In this regime, the relevant 
fields are the coordinates 4>i/ 4>n -{■!-, ■ ■ ■ , 4>n/4>n+i on CP N . They create BPS states of 
mass \m\. The field = 4>n + i/\/N plays the role of a Higgs field breaking the U(l) 
gauge symmetry (|2.13| ). Choosing the unitary gauge and using ( [4.10|) , we write 





\m\ \ 






V2tt 


W\) 



1/2 9 R 

g{\m\) 



where p is a fluctuating real scalar field. The constraint ( |2.12D , which is valid in 



the vacuum ( (4.10 ), shows that p is a composite operator creating a two-particle 



BPS/anti-BPS bound state of the elementary quanta. Though the attractive force 
between these quanta is of order 1/N at large N, the mixing between the N flavors 
will stabilize the bound state significantly, and the binding energy should be of order 
iV°. The mass of this bound state is a nice example of an observable which 
has a highly non-trivial perturbative expansion, as opposed to the cases studied in 
Section 3. We can straightforwardly calculate the leading large iV approximation for 
mb, by looking at the quadratic piece of the effective action deduced from (|4.7|) , see 
Appendix A. There is a mixing between f and the gauge multiplet.0 By inverting 
the matrix-valued propagator, we find a pole at p 2 = —ml such that 



In 



m 
A 



1 + Vl + u 2 1 8vr 2 8tt 2 „ n . 
arctan- = — — - , (4.12) 

3 A naive application of the standard results about the super-Higgs mechanism in four dimensions 
suggests that the Higgs and the gauge fields actually belong to the same supersymmetry multiplet 
and have the same mass. This is not correct in two dimensions, because a non-linear twisted 
superpotential is generated. 
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/ 4 m 2 


- 1 


U = t 






V m l 





where 

87T 2 /i 87T 2 , /i inN 

AJ~ lA|' AL^'-M' <4 ' 13) 

If we define perturbation theory in terms of the 't Hooft coupling constant A 2 with 
renormalization scale fi > |A|, ( |4.12j ) implies an expansion 

A 2 00 
8vr 



k=2 



with some /^-dependent coefficients Uk- Near the critical point m — > |A| or A^c — ► A^, 
Q4.12| ) implies 

M ~ l-A 2 /A 2 • (415) 

For the general picture of Section 4.1 to apply, one would need to prove that the series 
( |4.14j ) has a radius of convergence 

R» = K,c (4.16) 



or equivalently that 



A 



2(l-fc) 



k-^oo 87T 2 



(4.17) 



at least for a particular choice of /i. This would indeed imply that the perturbative 
series (|4.14|) is dominated by the terms with large k, or equivalently by the large 
Feynman graphs, near the singularity. Proving (|4.16|) , however, turns out to be 
particularly tricky. One can show rigorously that > for all /i, and that if 
Rfi — A MjC for some particular /z, then it is true for all larger values of /i as well. 
One can also show, using Picard theorem, that for /i = |A|, which means A^ iC = oo, 



the radius is actually finite, contradicting ( 4.16 ). To really understand what was 
going on, we performed a numerical analysis. For large enough values of /i, such 
that A 2 c /(87r 2 ) < 0.45, ( |4.17| ) is found to be satisfied, with a rapid and smooth 



convergence. For A 2 c /(87r 2 ) > 0.45, however, the behaviour of the coefficients 
changes drastically and ( |4.17| ) is apparently violated. 

The perturbative series for the mass of the bound state itself can be immediately 
deduced from ( [4.14j ) and fl4.13| ), and it has the same properties. At small coupling, 
we have 

m b ~ 2\m\ (l - + 0(AJ)) , (4.18) 

but near the critical points the high orders in perturbation theory dominate and we 
find 

m b oc >0. (4.19) 

u 
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It is important to realize that the asymptotic behaviours ( |4.15 ) or ( 4.19|) do not 



depend on /z, an obvious consequence of the renormalization group equations. If we 
introduce 5 = r — 1 ~ ln(m/A) as in ( |3.15| ), the asymptotics read 



u ~ 1/\S\ , m b ~ 2|A5| . (4.20) 

This equation has two important consequences. First it shows that the BPS/anti-BPS 
bound state becomes massless at the critical point in the leading large N approxima- 
tion. This is a nice result, because usually only BPS states can be proven to become 
massless at such singularities, by using exact BPS mass formulas as the one discussed 
in Section 2 or 3 (in our example the massless BPS states are the solitons interpolating 
between the two vacua that "collide" at the critical point). Second, after doing the 
rescaling ( f4.1|) for p = 1 (the correct value for a first class singularity), the mass m b 
is multiplied by N (in the same way as the central charge was mutiplied my N, see 
( P.16| )), and thus has a finite non-trivial limit in the double scaling ( |3.15|) at leading 
order, 

(iVm b ) scalcd = 2|A|/k + O(k ) . (4.21) 

What about the higher orders in 1/N? Can we trust the results obtained in the 
leading 1/N approximation? As we will show shortly, qualitatively, the answer is 
yes: the BPS/anti-BPS bound state does become massless, and Nm h does have a 
non-trivial finite limit in the scaling ( [3.15D . However, quantitatively, there are some 
important subtleties. The fact that (A^m b ) sca i e d is finite to all orders and has a non- 
trivial expansion in k (a result we will prove in the next subsection) implies that 
the 1/N corrections to m b must diverge at fixed N when 5 —>■ 0. The leading order 
equation (|4.20| ) is thus not to be trusted. The correct asymptotics is actually 

m b oc |A5 3/2 | . (4.22) 

One must not be confused and conclude that, in the exact theory, (|4 . 2 1|) is wrong. 



Equation ( |4.22| ) is valid when 5 — * at fixed N, while in the double scaling limit 
( p,15| ) we take N oo and 5 — > in a correlated way. The result ( |4 . 2 2| ) and the 
fact that (A^m b ) sca i e( i has a non-trivial expansion in k, far from being contradictory, 



actually complement each other. The non-trivial exponent in ( f4.22| ) is a consequence 



of the fact that the CFT at the critical point is non-trivial. This non-triviality is the 
cause of the divergences in the 1/N expansion. (A^m b ) sca i e d in turn picks up the most 
IR divergent terms in this expansion (see also for further discussion). 
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4.2.2 The double scaling limit 



Showing that the scaling ( |3.15| ) is fully consistent might look like a very difficult task, 
because it amounts to resumming the most divergent terms in the 1/N expansion to 
all orders and beyond. What makes it possible, and even easy, is the IR nature of 
the divergences. Not surprisingly, and as equation ( |4. 1| ) shows, this implies that the 
double scaling limit is also a low energy limit, and the path integral (|4.7| ) simplifies 
considerably in such a limit. The same property makes tractable the case of linear a 
models |23fl , and was first used in the context of non-linear a models in ||10|| . 



The starting point of the proof is the non-local effective action defining the path 
integral (fOP - By rescaling (J) — » (j)/y/N and introducing the field s = a + m and the 
functionals £ and ( discussed in Appendix A, it reads 

S eS = [d 2 xh a ^V a 4 > +(\s\ 2 -D)\4 > \ 2 + ^- *v + ^\n^- 

J V Z7T Z7T |A| 

+2N£[\s-m\ 2 -D,v a ] -2N( [-m + s,v a ] +fermions. (4.23) 



We will work thereafter with the bosonic fields only, the fermionic part of the action 
being unambiguously determined by supersymmetry. The rescaling of the space-time 
variables ( f4.1|) x a = Na a implies that a quantity of dimension D scale as N~ D . This 
means that the volume element d 2 x scales as N 2 , the partial derivatives d a and the 
fields s and v a scale as 1/N, and the fields D and v a p scale as 1/N 2 . Moreover, Q3.15D 



shows that 5 scales like 1/N. With those scalings, only a few terms in S e g survive 
when N — > oo. It is straightforward to check that those terms are at most linear in D 
and *v , and at most cubic in s. Terms containing derivatives cannot survive, because 
Lorentz invariance implies that derivatives must come in pair, and the dominant term 
with derivatives, J(ds) 2 , goes like iV -1 . These remarks imply that all the relevant 
terms in S e ff can be obtained from the potential ft4.8|) and from (|A.9|) . Adding up all 
the contributions, and using argm = argA + lm.5, we get 

S eS = [ cPx(v a ^V a <i>+(\s\ 2 -D)\ ( j)\ 2 + ^Re(5-s/A) 

J\A\ V 271 

_iN^v_ s _ + fermions \ + Q(ij N y (4 24) 



This formula is strictly valid only with a cut-off ~ A^°|A| that we have indicated 
on the integral sign, since we have been using a derivative expansion. The terms 
that potentially scale as iV 2 or iV cancel, which is a necessary condition for the 
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scaling ( 3.15 ) to be consistent. The terms cubic in s also cancel, as a consequence of 
super symmetry. Back to Minkowski space-time, and adding the fermions, ( |4.24|) can 
be written as 



S, 



cir 



d x 



d 4 9 &e 2V <S> -—Re I <W_<W, 
Air 



2A/ 



0(1/N). (4.25) 



It is useful at this point to introduce explicitly the scalings in N. This can be done 
in a manifestly supersymmetric way by defining 



x r 



Na>* , 8± = VN 6± 
which yield a new super field strenght 



6, 



N9- 



S 



ND + ^ x) D_^ x) V 



(4.26) 



(4.27) 



and an action 



S, 



cff 



d 2 a 

N\A\ 



d A Q $t e 2y $ 



Air 



Re / dQ_dQ+[ N5 S - 



2AJ 



+ 0(1/N). (4.28) 



The cut-off in the new space-time coordinates a is now of order N. Neglected terms, 
that all go to zero when iV — > oo, include for example the gauge and s fields kinetic 
terms that may be deduced from ( |4.23| ), 

N 



32tt|A| 5 



d 2 xd 4 9 SS 



1 

32tt|A| 2 A^ 



d'ad 4 B SS . 



(4.29) 



Let us now actually take the limit N — > oo and 5 —> in ( |4.28| ). Since the 



cut-off goes to infinity in this limit, one must renormalize the theory (|4.28|) in order 
to get a finite answer. This is the origin of the logarithmic correction to the naive 
scaling in ( |3.15| ). Only a one-loop renormalization of the linear term in the twisted 
superpotential is needed, as can be checked by integrating out the superfield $. One 
must add a counterterm S(\n(A / fi))/ (Air) to make the theory finite, where A = N\ A| 
is the cut-off for (J4.28[) and fi an arbitrary renormalization scale. This means that 



N5 is renormalized, with 



N5 = (N5) T + In 



A r 



(N8) x + In - — - + In A" . 



(4.30) 



/i fi 

where (N8) v = 1/k + InK is the renormalized quantity to be held fixed when the 
cut-off is removed. We thus recover the scaling ( p. 15 ), and this completes the proof. 

Several comments are here in order. First, it is important to understand the 
meaning of the "truncated" action (|4.28 ), with respect to the full action ( [4.23 ). To 
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do the ordinary 1/N expansion, one starts from ( |4.23| ) and expands around a saddle 
point, for example ( |4.9| ). An infinite number of vertices is then generated from ( |4.23| ), 
a vertex of order p contributing with a power of N 1 ~ p ^ 2 by the standard large N 
counting. The few terms that we have kept in ( |4.28| ) or ( |4.24| ), like the terms |s| 2 |0| 2 
or — D\<p\ 2 , correspond to the vertices producing the most IR divergent contributions 
near the critical points, which are the only one that survive in the scaling Q3.15| ). A 
second important comment is that the double scaled theory does not depend on a 
cut-off. It is a field theory consistent on all scales, defined by the action 



scaled 



d 2 a 



1 

47T 



Re / de_de_ 



— h In k + In — 

K jX 



s 



S 2 \ 



2k) 
(4.31) 

with the UV cut-off Ao taken to infinity and renormalized coupling constant k. Inter- 
estingly, the phenomenon of dimensional transmutation takes place, and the coupling 
k is actually replaced by a scale M is the quantum theory, 



— h In k 

K 



In — • 

M 



(4.32) 



The physics described by the action ( |4.31|) depends on the dimensionless ratio 

R = M/A . (4.33) 

By integrating out the superfield $, one can deduce the effective superpotential w by 
using (|2.371) , 



w(s) 



1 

47T 



s In 



eM 



s 

2A 



A 


u u 2 ' 




u In 


. ~~ 47T 


L eR 2 J 



(4.34) 



with s = Au. We recover, up to a global factor, the result obtained in Section 
3, equation ( p. 21 ). The double scaled theory has two vacua, obtained by solving 
dw/du = 0. When R = 1/e, the two vacua collide at u = 1, and we get a critical 
point, which is nothing but the original critical point used to define the double scaling 
limit. By expanding around u — 1, S — A + T, and by using the formulas of the 
Appendix, one can deduce the low energy effective action describing (|4.31 ) near the 
critical point, 



scaled, eff 



/ d 2 a 




[32n\A\ 2 j 


/|A| 



d 4 Q T^T Re 

47T 



d&_d& + (T\n(eR) + 



T*_ 
6A 2 



(4.35) 

Not surprisingly, we obtain a simple Landau- Ginzburg description of the A\ minimal 
M = 2 CFT. Note that the double scaled theory ( |4.31| ), however, differ from this 
simple Landau-Ginzburg description at high energies. 

Let us emphasize that the same qualitative phenomena are likely to occur in 
the gauge/string theory case at a first class singularity [§]. In particular, the string 
coupling should dissappear and be replaced by a mass scale. 
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4.3 Second class singularities 



For the sake of simplicity and conciseness, we will study the case of the simplest 
critical point only, corresponding to m = 2 in the notations of Section 3.3.2. Unlike 
the case of the first class singularity, we can integrate over all of the N + 1 superfields 
$i in (|4.5|), and we are left with the following path integral over S, 



V{a,v,D,v a )exv -(N + 1) J d 2 x In ^ + 



10 



2ir(N + l] 



* v 



det r (y - a - m) det r (y -a + rh) 



(4.36) 



N+l 
~ 2~~ 



det r (-V 2 + \<t + m| 2 - D - 2 V) det r (-V 2 + \a — m\ 2 — D — 2 JiFl 1 ^) 



where F± = y — (<7 ± m). The effective potential is 



D . fx |cr + m| z , | cr + m\ 



eff 



2 7 r' ln R + 



ln 



cr — m , la — ml 
+ J — ; L ln- 



57T 



a + m — D , a + m 



In 



2 -D 



8tt 



cr — m 



e/i 



In 



STT 



57T 



cr — m 



(4.37) 



and yields the iV — * oo vacuum expectation values 

(D)=0, \(a) 2 - m 2 \ = \A\ 2 . (4.38) 

The solution for a gives the enhangon that was described in the Figure 4 of Jl|. The 
critical point is obtained when the two disconnected components of the enhangon 
collide, at \m\ = |A| and (a) = 0. 

To study the double scaling limit, we use the same strategy as in Section 4.2.2. 
We focus on the bosonic part of the action. The rescaling of the space-time variables 
is x a = \^Na a (|4.1| ), showing that d 2 x scales as N, the partial derivative d a and the 
fields cr and v a scale as 1/yN, and the fields D and v a p scale as 1/N. As for the 
deviation from the critical point, 

2 

Tft 

6 = In , (4.39) 

equation ( ^.31[ ) shows that it scales as 1/N. Those scalings imply that the only terms 
that can survive are either the kinetic term for cr or the gauge field, that can be 
derived from ( |A.2C ) and ( |A.21 ), potential terms at most quartic in a and quadratic 



in D, that can be derived from (|4.37|) , and terms linear in v a p and at most quadratic 
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in a that can be derived from ( |A.9| ). Adding up all the contributions, we get the 
action 



S cS = (N + 1) / d 2 x I (\v aP v a p + d a a^d a a - D 2 ) 

y I a I o7T 777. \2 / 



|A| 

47T v ' 7 4vr 



Im(5 — cr 2 /m 2 ) + fermions 



+ 0(1/VN). (4.40) 



Introducing the scalings explicitly, 

x* = VN<t», # ± = Af 1/4 6 ± , # ± = Af 1/4 6 ± , S=VNZ, N5 = t, (4.41) 
( I4.4Q) reduces to, 



'scaled 



— / d ' z a 



1 



L4|A| : 



J d 4 Q SS + Re J de_de + (ts 



s 3 \ 



3A 2 



) 



(4.42) 



No remormalization is needed for the action ( |4.42|) , and thus the cut-off Ao = \/iV|A| 
can be removed. The double scaled theory for a second class singularity is thus given 
by a simple Landau-Ginzburg action, as was suggested in Section 3.3.2. It coincides 
at low energy with the first class double scaled theory ([4.35 ), but differs from it in 
the UV. 



5 J\f = 1 supersymmetric models 

The purpose of the following brief Section is to emphasize the fact that the results 
obtained so far do not depend on supersymmetry. We focused on a supersymmetric 
model because our main goal was to make a comparison with the four dimensional 
supersymmetric gauge theories studied in |J 0. In fact, we believe that constructions 
of non-supersymmetric four dimensional non-critical strings could be made by using 
non-supersymmetric gauge theories with Higgs fields and adjusting the parameters in 
the Higgs potential to approach a critical point. 

The M = 2 supersymmetric model we have studied in details in this paper was 
based on the classical potential ( 2.21 ) with the constraints ( |2.12| , |2.13| ) on the complex 



fields 4>i. It is natural to suspect that a very similar, but only J\f = 1 supersymmetric, 
model could be constructed for which the classical potential is 



j N+l 

v cl = -Y,(° + m) 2 <f>h (5.1) 
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with real fields fa, a and mass parameters m,, and constraint 

N+l 

E*f = £- («) 

i=l y 

Such a model indeed exists. By introducing Majorana spinors ipi which are in the 
same supermultiplet as the fas and a supermultiplet (//, .D) of Lagrange multipliers, 
its lagrangian reads 

N+l 



+f (E *E**- < 5 ' 3 > 

i=l ^ i=l 

This is a £? non-linear cr model with M — 1 supersymmetry. The mass terms come 
from a superpotential 

W(^) = -J2^l (5.4) 

8=1 

The lagrangian ( |5.3|) is very similar to ( |2.20|) , and the large A" limit can be studied 



with the same methods [24]. A particularly interesting aspect of the model ( |5.3j ) 
is that the number of vacua changes when the mass parameters are varied, while 
the Witten index tr(-l) F = 1 + (-1)* is, of course, constant. In the CP N model, 
holomorphicity was preventing such drastic changes in the space of vacua. At large 
N, one can show that ( [2.35D is replaced by 

N+l 

n^ + m,) 2 = AW), (5.5) 

i=l 

where all the parameters are now real. In the A — > weak coupling limit, we have 
2 (A" + 1) vacua coming 'm N + 1 inequivalent pairs, 

(o-)w.c. ^ -m , (4>j) w . c . ^ ±h^i , (5.6) 

while at strong coupling we are left with only two vacua 

(<7> BX . ~±A. (5.7) 

The particular values of the mass parameters for which the number of vacua changes 
correspond to critical points where the large A" expansion breaks down, and where 
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double scaling limits can be defined. The double scaled theories are typically simple 
Landau- Ginzburg theories [2^]. The scalings always look like ( |3.15| ), with logarithmic 



correction, because M = 1 Landau- Ginzburg models need to be renormalized by 
normal ordering the superpotential. 

Instead of considering superpotential-induced mass terms, one can also use the 
isometries of the target space, as was done for the CP N model. In that case, the 
M = 1 supersymmetric lagrangian takes the form 

N+l ^ N+l 



k,l=l i=l 

i=l " i=l 

where is an antisymmetric mass matrix. Again, critical points can be found, and 
double scaling limits defined pifl . An interesting aspect of the model ( p.8| ) is that 
supersymmetry can be spontaneously broken when N is odd. 

In both models ( |5.3| ) and (p.8|), the masses induce a quadradic term hij(j>i4>j/2, 
where the tensor "magnetic field" h = mm T is expressed in term of a symmetric 
or antisymmetric mass matrix respectively. One cannot find a supersymmetric the- 
ory with arbitrary h, but the corresponding non-supersymmetric model can also be 
studied, and again critical points are found and double scaling limits can be defined 

m 



6 Prospects 

The main goal of the present paper was to improve one's understanding of the results 
obtained in |l|, f§. We hope that our analysis has convinced the reader that the 
gauge theories double scaling limits are likely to yield well-defined four dimensional 
non-critical string theories, as conjectured in 0. A natural avenue for future work 
is to try to generalize the cases studied in @. There is a variety of critical points 
appearing on the moduli space of supersymmetric gauge theories, and a large class 
of string theories can certainly be generated. Remarkably, for all those theories, the 
dimensions of the world sheet couplings as well as the space-time central charge as a 
function of these couplings can be calculated exactly. It would be interesting to study 
in details the structure of the formulas so obtained. One may hope, from experience 
with the c < 1 matrix models where the KdV hierarchy plays a prominent role 
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22 1 , that a general mathematical structure could emerge. Unravelling this structure 
might eventually help in understanding our string theories from a more conventional, 
'continuous' point of view. 

Another fascinating possible direction of research is based on the fact, explained 
in the Introduction, that our two dimensional models admit brane constructions. A 
crucial feature is that the large A" limit of the a models correspond to a large number 
of branes, as in the case of gauge theories PS[ . One might then expect to be able 
to find a description involving quantum gravity. The startling point is that, in sharp 
contrast with the gauge theory case, the large A" limits of our models are exactly 
solvable. 
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A Formulas for determinants 

We consider the renormalized euclidean determinants 

det r (-V 2 + h)= exp(2£[Ma]) , det r (y - /) = exp(2 C[j> a ]) , (A.l) 
where the functionals £ and ( are defined by the equations 

= itrln(-V 2 + h)-£ [h,v a ], (A.2) 

C[/, v a ] = ~ tr ln(y - /) - Co\f, v a ] , (A.3) 



with local counterterms £o and The covariant derivative is V a = d a +iv a , his & real 
and positive field, and / is a complex field with associated matrix / = Re / — ia 3 Im /. 
The local counterterms depend on the regularization and renormalization schemes. 
Gauge invariant and supersymmetric results can be obtained by using a Pauli-Villars 
regularization with cut-off Ao, renormalization scale \x and counterterms 

£o = f fd\h, Co = ^ln^ [d*x\f\ 2 . (A.4) 

47r n J An fi J 
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Dimensional regularization D = 2 — e can also be used, with the definition 

± hl h = /j^_L_ = _L_J_ hl £_:L, (A .5) 
2f /i j (2n) D p 2 + /i 2 2vre 4tt 4tt 4tt v ; 

but one must add a finite local counterterm — J d 2 x(lm f) 2 / (4tt) to the fermionic 
functional for supersymmetry to be preserved (such a term is generated due to the 
unusual properties of the tjja^iplva f vertex in dimensional regularization). 

A.l Special cases 

For v a pure gauge, h constant and positive, and / constant, we have 

£[h = est > 0,v a = d aX ] = J d 2 xV h {h) , C[/ = cst,t; a = d aX ] = J d 2 xV { (f) , (A.6) 

with the potentials 

KW = -^ln^, W) = VUI/| 2 ). (A.7) 

If / = 0, we have, for any v a p going to zero fast enough at infinity , 

If 1 

C[/ = 0, v a ] = — / d 2 x v a(3 — v af3 . (A.8) 
If Im// Re f is constant, the term linear in v a p in ( can also be exactly calculated, 
({f, vjiincar in fi = ^ J d 2 x e aP v a p Im In / • (A. 9) 

A. 2 General case 

In general, one writes 

h = m 2 + ip, f = M + <f>, (A.10) 

with m real and M complex, and one expands the functionals in powers of the fields 

<f, (p and v a , 

oo oo 

t[h,v a ] = ^2£ n [<p,v a ;m 2 ], ([f,v a ] = ^ Cn[& v a ; M] . (A.ll) 

n=0 n=0 

The functionals 

& = -Lin 4 f<P*<P, (A.12) 
87T m z J 

Ci = ^ j d 2 xe af3 v af3 Im\nM + ^ hi ^ J ' d 2 x Re(M» , (A. 13) 
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are linear in the fields, s 2 and (2 are quadratic, etc. . . The quadradic pieces are most 
easily expressed by introducing the Fourier transforms of the fields, 



= Jd 2 xe ikx v a (x) , etc, 



and the one-loop integral 

S{k;m 2 ) 



d 2 p 



8tt 2 J (p 2 + m 2 )((p+k) 2 + m 2 ) 



(A.14) 



(A.15) 



This integral can be evaluated in different regimes (euclidean, and below or above the 
pair production threshold), by using Feynman's ie prescription when necessary, 



S(k; m 2 



ln^ 1 + 4m2/p2 + 1 , for ^>0, 

4vr p z y/l + Am 2 /p 2 a/1 + Am 2 /p 2 - 1 



2ir p 2 ^J-l - Am 2 /p 2 



arctan 



1 



a/— 1 — Am 2 /p 2 



for < -p 2 < Am 2 , 



1 1 . 1 + a/1 + 4m 2 /p 2 . \ 2 A 2 

7: — , In , — in , tor — zr > Am A . 

4?r p 2 - s /l + 4m 2 /p 2 ^ 1 - a/1 + 4m 2 /p 2 / 

(A. 16) 

Introducing 0i(x) = Re0(x), 02 (a;) = Im0(x), Mi = ReM and M 2 = ImM, we have 



1 f d 2 k 



c 2 = 



2 J (27if 



rf 2 A; 



(A; 2 + 4m 2 )^(A;;m 2 



1 \ /fc a fc/3 



2tt/V jfe 2 



+S(k;m 2 )ip(-k)(p(k) 



, (A.17) 



2 J (2tt) ; 



(4|M| 2 S(A;; |M| 2 ) + i-) - <w)^(-^W 

+ ((A; 2 + AM 2 )S(k; \M\ 2 ) + ±- In ^^(-fc)^*) 

+ ((A; 2 + 4M 2 )S(A;; |M| 2 ) + i- In jjL^-fc)^) 
+AS(k; |M| 2 )e a/3 A; /3 ^(-A;)(M 1 2 (A;) - M 2 <Mfc)) 



+8M!M 2> S(A;; |M| 2 ) 0i(-fc)0 2 (fc) 



(A.18) 



The low energy expansion up to two derivative terms is obtained by using 

k 2 



S(k;m 2 ) 



(A.19) 
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and it reads 



I.e. 



487rm 2 



487t|M| 2 



1 



1 




d 2 x -v a pv a p + -^—^d a (pd a p - 3y? 2 , 
/ d 2 x -v af3 v af3 + Qi lm(M*(f))e al3 v al 3 



3 d a $d a <i> 



(A.20) 



2 



Re{Md a <p) Re(Md a (j)) + l2\M<j)\ 2 In 



12 (Re(M*0)) 2 



(A.21) 
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